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ABSTRACT 

We apply a von Zeipel gravity darkening model to corotating binaries to obtain a simple, analytical 
expression for the emergent radiative flux from a tidally distorted primary orbiting a point-mass 
secondary. We adopt a simple Roche model to determine the envelope structure of the primary, 
assumed massive and centrally condensed, and use the results to calculate the flux. As for single 
rotating stars, gravity darkening reduces the flux along the stellar equator of the primary, but, unlike 
for rotating stars, we find that gravity brightening enhances the flux in a region around the stellar 
poles. We identify a critical limiting separation beyond which hydrostatic equilibrium no longer is 
possible, whereby the flux vanishes at the point on the stellar equator of the primary facing the 
companion. For equal-mass binaries, the total luminosity is reduced by about 13% when this limiting 
separation is reached. 
Subject headings: binaries: close - stars: rotation 



1. INTRODUCTION 



Early studies (von Zeipel 1924a|b|c Ghandrasekhar 



1933 1 showed that on stellar surfaces the radiative hux is 
proportion al to the effective gravitational force (see, e.g., 
[Kippcnhahn fc Weigert_1990j Tassoul_2000^ for reviews). 
In rotating stars, this means that regions close to the pole 
are brighter (and have a higher effective temperature) 
than regions close to the equator, an effect that is often 
referred to as gravity darkening. In fact, for stars rotat- 
ing at the break-up s peed, the radiative flux vanis hes at 
the equator (see, e.g., Baumgarte fc Shapiro||19 99i here- 
after Paper l) . Gravity darkening plays an imiportant role 



in the classification of stars (e.g. Maeder & Peytremann 



19701, and may even provide a means of estimat ing stel 
far mass es independently of binary companions ( Zhao et 



al 



hile the simplest gravity darkening models make sev- 
eral assumptions that restrict their applicability to very 
massive and supermassive stars (as we will review in Sec- 
tion [2] below) , the y continue to be useful in a number of 
different contexts. |Granmer fc Owocki (1995), for exam- 
ple, adopt similar models to compute radiative fluxes and 
driving forces for winds from rapidly rotating B stars. In 
Paper I, two of us adopted a simple Roche model to de- 
termine the envelope structure of and emergent radiation 
flux from a rotating, supermassive star (SMS) and em- 
ployed the results to calculate the evolutionary timescale 
for these objects. 

An exciting recent development provides a new moti- 
vation for studying gravity darkening. With optical or 
infrared interferometric arrays (including the PTI, NPOI 
and GHARA arrays), it has become possible to obtain 
resolved in terferome tric images of individual stars (e.g. 
van Belle et al. 200ll [Peterson et al.T2006l [Aufdenberg et| 



al. 2006). Observations of rotating stars, which indeed 
show oblate shapes with brighter poles and darker equa- 
tors, are analyzed using gravity darkening models. These 
models can be used to relate different stellar features to 
each other, for example the ratio between the polar and 



equatorial fluxes to the star's angular velocity. While 
some deviations between observations and models point 
to the fact that some of the assumptions made in the 
models may be too restrictive (for example, the observed 
stars may have c omplicated surface lay ers or differential 
rotati on; see, e.g., Mq nnier et al. ( 2007 ) as well as Section 
7 of Zhao et al. ( 2009p for discussions), the basic features 
are represented reasonably well. Moreover, the quality of 
fits can be improved by introducing additional free pa- 
rameters into the models (for other model im prov ements 
see, e .g., [Espinosa Lara fc Rieutord (2011) and Claret 
( [2012| ). 

In addition to observations of rotating stars, resolved 
interfero metric images of c lose binary stars have become 
available ( Zhao et al.|2008 |. While these observations do 
not yet have sutticieiit resolution to distinguish detailed 
models, they do motivate a study of gravity darkening in 
binaries. 

In this short paper we point out that, under the as- 
sumptions listed in Section [2] the simple gravity dark- 
ening models for rotating stars can be generalized very 
easily to describe corotating stars in a binary system. 
As a consequence, we obtain very simple analytical ex- 
pressions for the envelope structure, flux and the total 
luminosity for a massive or supermassive primary orbit- 
ing a point-mass secondary as a function of the binary 
separation and mass ratio. While the effects of rotation 
always decrease the surface flux (compared to a nonro- 
tating star in isolation), the presence of a binary com- 
panion decreases the flux in some regions of the star, but 
increases the flux in a region around the pole. In a bi- 
nary, we therefore encounter both gravity darkening and 
gravity brightening. We also identify a critical limiting 
separation below which hydrostatic equilibrium cannot 
exist and at which the surface flux from the point on 
the stellar equator facing the binary companion vanishes. 
For an equal-mass binary, the total luminosity is reduced 
by about 13% once this separation has been reached. 

Elements of this problem have been known for a long 
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time (e.g. von Zeipel 1924c |Osaki 1965) and the model 
we construct presumably will be superseded by more de- 
tailed numerical treatments. We nevertheless hope that 
our simple analysis and analytical expressions may prove 
useful, for instance, in providing preliminary interpreta- 
tions of future resolved images of close binary stars and 
calibrating more detailed numerical models. 

Our paper is organized as follows. In Section [2] we 
review our basic assumptions. In SectionlSlwe review the 
Roche approximation and its application to binaries to 
determine the envelope structure and locate the surface 
of corotating stars in a binary. In Section [4] we compute 
the radiative flux from these stars, and find both gravity 
darkening and brightening. We integrate this flux to find 
the total luminosity in Section [5] and conclude with a 
brief discussion in Section [H 

2. BASIC ASSUMPTIONS 

Our analysis of a binary relies on several explicit 
assumptions. Some of these assumptions we already 
adopted in Paper I, where we treated uniformly rotating, 
highly massive stars and SMSs in isolation. To determine 
the equilibrium structure of the envelope, we assume that 
the primary is 

1. governed by Newtonian gravitation, 

2. in synchronous orbit about a point-mass secondary, 
and that its envelope 

3. is described by the Roche model, 

4. is characterized by a polytropic equation of state, 
and 

5. interacts with the secondary via a potential that 
can be truncated beyond the quadrupolar tidal 
term. 

To determine the emergent radiative flux from the pri- 
mary, we further assume that its envelope is 

6. dominated by thermal radiation pressure, 

7. fully convective, and 

8. characterized by a constant Rosseland mean opac- 
ity (e.g., electron scattering). 

We assume that gravitational fields are sufficiently 
weak so that we can apply Newtonian gravity. This as- 
sumption clearly holds for normal stars. SMSs of great- 
est astrophysical interest have masses and radii that sat- 
isfy R/M > 400 so that this assumption certainly holds. 
Relativistic corrections are important for the stability of 
SMSs, but can be neglected in the analysis of the equi- 
librium state. 

The assumption that the binary is in corotation follows 
from the observation that most short period binaries con- 
taining massive stars orb it in circular, synchronous orbits 
( Vanbeveren et al.|1998 ) . The combination of convection 
and magnetic fields are likely to generate an effective tur- 
bulent viscosity, which dampens nonsynchronous motion 
and brings the binary into corotation. 



The assumption of a point-mass companion is adopted 
for simplicity. If the corotating companion is also a mas- 
sive star, the point-mass can be replaced by a finite star 
whose envelope can be treated identically to the enve- 
lope of the primary. The results for the companion's 
envelope structure and emergent ffux would depend the 
same way on its mass and radius as we find for the pri- 
mary. We truncate the potential of the secondary after 
the quadrupolar tidal term, which captures the leading- 
order, dominant effects of the secondary on the primary. 

For large masses, the ratio between radiation pressure, 
Pr, and gas pressure, Pg, satisfies 



Pr \Mq 



-1/2 



(1) 



see , e.g., eqs. (17.2.8) and (17.3.5) in Shapiro fc Teukol^ 



sky 1983); here the coefficient has been evaluated tor a 
composition of pure ionized hydrogen. For SMSs with 
M > 1O^M0, we can therefore neglect the pressure con- 
tributions of the plasma in determining the equilibrium 
profile, even though the plasma ma y be important for de- 
termining the stability of the star ( Zel'dovich & Novikov 
19711 IShapiro & Teukolsky''1983|). A simple proof that 



very massive stars or SMSs are convecti ve in this limit 
has been given by Loeb & Rasio ( 1994[ )p] This result 
implies that the photon entropy per baryon. 



4 aT^ 
3 ub 



(2) 



is constant throughout the star, and so therefore is f3p w 
8{sr/k)~^. Here a is the radiation density constant, T 
is the temperature, n^ is the baryon density, and k is 
Boltzmann's constant. As a consequence, the equation 
of state of a very massive star or SMS is that of an n = 3 
polytrope: 



P = Kp^'\ K = 



k 



5{l + (3pf 
a p% 



-,1/3 



= const. 



(3) 

where m is the atomic mass u nit and /i is the mean molec- 
ular weight (cf. Clayton 1983, eq. (2-289); note that Clay- 
ton adopts a ditterent definition of /3p, which is related 
to ours by ^ciayton = Pp / {I + Pp)) . In the high tem- 
perature, low density, strongly ionized plasma of a very 
massive or SMS, Thomson scattering off free electrons is 
the dominant source of opacity. This opacity is indepen- 
dent of density and justifies our assumption about the 
Rosseland mean opacity. 

In applications to SMSs our analysis neglects electron- 
positron pairs and Klein-Nishina corrections to the 
electron-scatteri ng opacitv. which is valid for Af > 
IO^Mq (see, e.g. [Fuller et"ar]|l986| ). 

3. THE ROCHE MODEL FOR A BINARY 

We begin with the equation of hydrostatic equilibrium 
satisfied by the primary. 



VP 



^V($p + $c + $r). 



(4) 



1 In fully convective stars energy transport in the stellar interior 
is dominated by convection. However, there is still an outgoing flux 
of radiation from the stellar surface. The later is what we compute 
in Section 4 below. 
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Fig. 1. — Sketch of the coordinate system used in our calculation. 
The primary of mass Mp is represented by the shaded configuration 
on the right, while the companion of mass Mc is represented by 
the dot on the left. 

Here the right-hand side describes an effective gravita- 
tional force, which is derived from the (interior) gravita- 
tional potential ^p of the primary, the (exterior) gravi- 
tational potential of the companion ^c, and a centrifugal 
potential ^.^ arising from the (synchronous) rotation of 
the primary. 

Stars with soft equations of state are centrally con- 
densed, i.e. most of the mass is concentrated in a high- 
density core that is surrounded by an extended low- 
density envelope. For an n ~ 3 polytrope, for ex- 
ample, the ratio between central and mean density is 
Pc/p = 54.2. The gravitational force in the envelope is 
therefore dominated by the massive core, and it is thus 
legitimate to neglect the self-gravity of the envelope. In 
the envelope, we may therefore approximate the Newto- 
nian potential of the primary star <I>p as 



*„ = -- 



Ah 



(5) 



where Mp is the mass of the primary and r the distance 
from the primary's center (here we adopt gravitational 
units by setting G = 1). 

In accord with the Roche model, the companion may 
be treated as a point mass M^, located at a distance s 
from the center of the primary (see Fig.ll]). The potential 
at a point r in the primary is then given by 



$c = -- 



Mr. 



(6) 



where d is the distance from Mc- We now expand d about 
s > r, 

where Pe{cos9') is the Legendre polynomial of order £, 
and where 6' is the angle between r and the line con- 
necting the center of the primary to the point mass com- 
panion. The first term ^ = is a constant that can be 
ignored. We will see that the £ = 1 term will cancel out 
later when we consider the rotational contribution <I>r, 
but we retain it for now. The first relevant term is the 
quadrupolar tidal term £ — 2. Truncating the expansion 
(l7| after this term we approximate 



$. 



Mrr 



cos 9' 



M^ 
2s3 



(Scos^ei'-l) 



(8) 



We now introduce a coordinate system as shown in Fig.fTl 
so that the orbital plane is in the x — y plane, with the 
center of the primary at the origin and the companion at 
X = —s and y = z — 0. In terms of spherical polar co- 
ordinates we now express cos 9' = —x/r = — cos0sin0, 
and hence 



$. = 



MrX 



^t(3cos 



)Sin 



-i). 



(9) 



s2 2 s 

Finally, the rotational potential ^^ in Q arises from 
the rotation of the primary about the system's center 
of mass. Assuming corotation, so that the star appears 
static in the rotating frame of the binary, we may write 
this term as 



-^r = -^^^{{x-xcMf+y'^), 



where 



Mr 



a;cM 



M^ + Mr 



(10) 



(11) 



is the location of the center of mass. Using (11) as well 
as Kepler's law 

Mp + M, 



n' 



we can write ( 10 1 as 



$, 



_ 1 2 McX 

— "2 CM -^ 



1 Mr> + Mc 



r^ sin^ ( 



(12) 



(13) 



The first term on the right-hand side is a constant term 
that can be ignored. The second term will exactly cancel 
the first term in Q when the two potentials are added 
in Q ; wcmay therefore discard the linear terms in both 
([9| and (13). Our potentials $£ and ^r then reduce to 



^r + ^r 



-^4(3cos2 0sin2 0- 
2 s"^ 



1) 



IMp 
2 



Mr 



2 • 2 

r sm I 



Integrating eq. ^ yields the Euler integral 

h + ^p + ^c + '^T = H, 
where H is a constant of integration and where 



h = 



dP , ,P 

— = (n+l)- 
P P 



(14) 



(15) 



(16) 



is the enthalpy per unit mass. Evaluating ( |15[ ) for infinite 
binary separation s — >■ oo (where $c = ^r — 0) at the 
stellar surface (where /i = 0) we find 



H 



Mp 



(17) 



where Rq is the stellar radius for the nonrotating (spher- 
ical) star in isolation. 

It is consistent with the Roche approximation to as- 
sume that the core remains unaffected by the presence 
of the binary companion. Given that $c and $r vanish 
at the center of the star, and given that we may assume 
h and $p to remain unperturbed there, we may also as- 
sume that H is independent of the binary separation, 
so that we can always express it in terms of Rq as in 
eq. (17). Similar arguments have been made for isolated 
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Fig. 2. — The dimensionless stellar radius z as a function of cosS 
for an equal-mass binary (q = 1) at different values of the binary 
separation a. For each binary separation the top line represents 
results for (j> = n (i.e. in the direction toward the binary companion) 
while the bottom line represents results for <f> = 7r/2. The two lines 
connect at the pole (cos 9 = 1). Note that for a = criim, w e ha ve 
z = 3/2 for the point facing the binary companion (see eq. ||26|l). 



rotating stars (e.g . Zel'dovich & Novikov 1971 Shapiro 



& Teukolsky 1983), where they have been conlirmeci by 
numerical simulations (e.g. Papaloizou fc Whelan||1973 ). 
The surface of the star satishes h = and can be de- 
termined from cq. (|15|) in the form 



(19) 



$p + $c + $r - if = 0. (18) 

Introducing dimensional parameters for the mass ratio 



the binary separation 



a 



s 



(20) 



and the distance from the stehar center to the surface 

r 

Kq 



(21) 



we can bring (18) into the form of a cubic equation for 

Csz^ - z + 1 = 0, (22) 

where we have abbreviated 
1 



C. 



2cr3 



(q{3 cos^ (f> sin^ 6* - 1) + (1 + g) sin^ o) . (23) 

Given binary parameters q and cr, this cubic equation 
can be solved analytically to find z as a function of the 
coordinates 9 and (j) (s ee A ppendix). 

Note that the cubic (22) takes exactly the same form 
as that for a single rotatmg star (see, eq. (23) in Paper I), 
except that the coefficient C3 is now different jj Instead 
of depending on 6 only, as for axisymmetric, rotating 

^ Our coefficient C3 reduces to the corresponding coefficient for 
rotating stars when 5 = 0. 



stars, it now depends on both 9 and 0. Moreover, for ro- 
tating stars the corresponding coefficient is nonnegative, 
resulting in values for z that are always greater than or 
equal to unity. Here, however, C3 can be positive or neg- 
ative. At the poles, where sin0 = 0, for example, we 
have C3 < (for finite a), resulting in z < 1. This is 
consistent with the fact that the tidal field of the com- 
panion leads to the squeezing of the primary along its 
poles and an elongation along its equator. 

Hydrostatic equilibrium is only possible if, at the sur- 
face of the star, where the density vanishes, the pressure 
increases towards the interior of the star. Our sequences 
of hydrostatic equilibria therefore terminate when, at 
the point facing the companion, i.e. at 9 = 7r/2 and 
(j) = TT, the right-hand side of Q vanishesP] Evaluat- 
ing the right-hand side of Q at that point, and setting 
it to zero, yields 



For 6*: 



(l + 3g)— = 1. 

tt/2 and cj) — tt we also have 
l + 3q 



Cs 



2ct3 



(24) 



(25) 



z, at 



Inserting (24) together with (25|) into the cubic (|22|) we 

find that the limiting value of 

is 



^lii 



7r/2 and <p — t:, 
(26) 



Interestingly, this is the same value obtained for the equa- 
tor of single stars rotating at the break-up limit (see, e.g., 
eq. (10) in Paper I). Inserting (26 ) back into ( p4| we now 
obtain the limiting binary separation crum for which, un- 
der our assumptions, sequences of hydrostatic equilibria 
end, 

^iim = ^(l + 3g)i/3. (27) 

For a mass ratio of g = 1, for example, we have CTiim ~ 
2.38. Given that a > 2.4 for all equilibrium models, our 
truncation of the interaction potential beyond the tidal 
term is justified in a first approximation. The corre- 
sponding Roche limit for a homogeneous, incompressible 
{n = 0) star, allowing for departure of the angular veloc- 
ity from the Keplerian value due to the ellipsoidal s hape 
of the primary, yields aiim = 2.713 ( Lai et aL]|1993 ). 

In Fig. [2] we show results for z for an equal- rnass binary 
for different values of the binary separation. In particu- 
lar, these results confirm the limiting value ( 26 ) for the 
end point of the sequence at crum. 

4. GRAVITY DARKENING AND BRIGHTENING 

In the diffusion approximation, the radiation flirx is 
given by 

where U is the energy density of the radiation, 

U ^aT^ ^ 3P, (29) 

•^ Having truncated the expansion ([tI after the tidal term, this 
point is equivalent to the point poinding away from the binary 
companion, 9 = 7r/2 and = 0. 
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Fig. 3. — The flux -F/-FEdd ^^ ^ function of function of cos 6 
for an equal-mass binary (q = 1) at different values of the binary 
separation a. For each binary separation the bottom line represents 
results for (j> = n (i.e. in the direction toward the binary companion) 
while the top line represents results for <j) = tt/2. The two lines 
connect at the pole (cos 6 = 1). Note that in a region around the 
pole the flux exceeds the corresponding Eddington flux, and that 
the flux vanishes at the point facing the binary companion when 






and where k is the opacity (which we assume to be dom- 
inated by electron scattering, k = Kos)- We have also as- 
sumed that the pressure is dominated by radiation pres- 
sure 

Pk.P^^ -aT^. (30) 

Inserting equations ([28|_and ( |29| ) into the equation of 
hydrostatic equilibrium (H) yields 



kF = V($p + $c 



<^r 



(31) 



In polar coordinates in an orthonormal basis, the mag- 
nitude of the flux is 



F = (F| 



F2 



Flf"- 



(32) 



Evaluating the gradients of the potentials <I>j,, $c and ^r 
wc find 



F 



F 



Edd 



((> 



^tZ^(3cos^ (f)sm^ I 



■^z^{3 cos^ (p COS 9 sin ( 



1 



9 3-2 



1 



q r, 

:—z sin cos ( 



( — z'^(3cos0sin0sin(?) J > 



1/2 



where 



^Edd — 



M, 



nr 



(33) 

(34) 



is the Eddington flux from a spherical star of radius r. 

We graph the flux F for an equal-mass binary for dif- 
ferent values of the binary separation a = s/Rq in Fig. [3] 
In a region around the poles the flux is greater than trie 
Eddington flux. The reason for this "gravity brighten- 
ing" effect is that, at the poles, the tidal forces caused by 



Fig. 4. — Surface images of the primary in an equal-mass binary 
at three different binary separations cr = 5 (top panel), cr = 3 
(middle panel) and a = g-iim (bottom panel). The shape of the 
prim ary is given by eq. | |22| , while the color coding represents the 
flux ( |33[ l. Yellow indicates the largest flux, red a smaller flux, and 
black a vanishing flux. The dot represents the point mass that 
models the companion. 

the companion lead to an increase of the effective gravi- 
tational force, which in turn leads to an increase in the 
flux. The opposite is true at the points on the stellar 
equator either pointing directly toward the companion 
or directly away from the companion. At these points, 
the tidal forces lead to a reduction in the effective gravi- 
tational force. On the equator, stellar rotation also leads 
to a reduction of the effective gravitational force. At the 
two points facing towards or away from the companion 
(i.e. ioi (p — TT or (p — 0) , the two effects act together to 
result in the greatest reduction in the brightness, i.e. the 
strongest gravity darkening effect. For the limiting bi- 
nary separation a — crum, the flux completely vanishes 
at those points. At the two points on the equator with 
(j) = tt/2 and = 37r/2, on the other hand, the two effects 
counteract, leading to a smaller reduction in the flux. 

In Fig. HI we also show surface images of the primary 
in an equal-mass binary at three different binary sepa- 
rations. The shape of the primary is given by eq. (22 1 
and reflects the tidal deformation, while the color coo- 
ing indicates the local radiative surface flux ( 33 ) . Yellow 



indicates the largest flux, red a smaller flux, and black a 
vanishing flux. 

5. LUMINOSITY 

We find the total luminosity of the star by integrating 
the flux ( 33 1 over the stellar surface 



F-dA 



(b FdA. 



(35) 
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Fig. 5. — The luminosity L/L^^^ as a function of binary separa- 
tion a for different values of the mass ratio q. For each valu e of q, 
the sequences start at cr = 10 and end at CTuni given by cq. (|27|). 



The surface element can be written as 



dA= il 



z^ sin I 



1/2 



^dn, 



^(36) 
where dfl = sm9d6d(j). We then insert the flux |33| and 
evaluate the integration numerically to find the star's 
luminosity L. It is convenient to divide the result by the 
Eddington luminosity 



L 



AttM 



Edd 



(37) 



and express the results in terms of the dimensionless ratio 
L/L-Edd- In Fig. p^ we show results for the luminosity as 
a function of the oinary separation a for different values 
of the mass ratio q. All sequences shown start at cr = 10 
and end at the limiting binary separation aum given by 
eq. (27). For an equal-mass binary, for example, the 
luminosity is reduced by about 12.7% when the binary 
reaches criii„. We also verified that in the limit g — ?> 0, 
the luminosity for a — a\i„^ is reduced by about 36%, 
which is the value found for single stars rotating at the 
break-up limit (see, e.g.. Paper I). 

6. DISCUSSION 

We apply gravity darkening models to corotating bi- 
nary stars and obtain simple and analytical expressions 
for the surface flux of tidally distorted stars. The tidal 



interaction leads to both gravity darkening (along the 
equator) and gravity brightening (in regions around the 
poles). We identify a critical separation at which, within 
the Roche model, sequences of hydrostatic equilibrium 
end, and at which the radiative flux at the point on the 
equator that faces the binary companion vanishes. At 
this critical separation, the total luminosity in an equal- 
mass binary is about 13% less than for the corresponding 
nonrotating star in isolation. 

Simple and analytical models for the flux and lumi- 
nosity from binary stars, even if they are approximate, 
are useful in many ways. In particular, we hope that 
they will provide useful models for comparison with fu- 
ture resolved interferometric images of close binary stars. 
For single rotating stars, von Zeipel gravity darkening 
models capture the basic features of interferometric im- 
ages reasonably well, but they also show some devia- 
tions. Presumably, these deviations are caused by the 
fact that some of the assumptions do not apply to the 
observed stars. In particular, many main sequence stars 
have complicated surface layers or rota te differentially 
rather than uniformly. As discussed by |Monnier et al.| 
( 2007 1 (see also Zhao et al. 2009 ) , the agreement between 
model and observatio ns can be improved by introducing 



new free parameters. Monnier et al. (2007) and others 



have used "/3-free" models, in which the ettective temper 
ature is taken to be Toff ex g^f^, where g^e is the effective 
gravitational force on the right-hand side of eq. ^ . In 
the "standard" model, (3 = 1/4, but allowing /3 to be 
a free parameter allows for improved fits to the obser- 
vations. Other improvements t o the simple von Zeipel 
models hav e be en proposed b y [Espinosa Lara fc Rieu-| 
tord| ( |2011[ ) and |Claret| ( |2012[ ). We expect that similar 
generalizations may improve future fits between binary 
gravity darkening models and observations as well. 
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APPENDIX 
SOLUTION FOR THE STELLAR SURFACE 

In this brief appendix we present the solution to the cubic equation (|22[ ) , yield ing the stellar surface z. The general 
solution t o a cubic equation can be found, for example, in Press et aT (2007). Applying their prescription to our 
equation (22), we see that the form of the solution z depends on the coefficient C3. 

By combining equations (23) and (27) we first observe that we always have C3 < 4/27. The form of the solution 
then depends only on the sign of Ca.Ti C3 is positive, then the cubic has three real roots and we pick the one that 
yields z — 1 when C3 = 0, given by 



arccos 



: cos 



(^27^74) - 2^ 



(Ca > 0). 



(Al) 



GRAVITY DARKENING AND BRIGHTENING IN BINARIES 



For a single, rotating star, we always have < C3 < 4/27 (see, e.g.. Paper I), so that the solution can always be 
written in this form. For a binary, however, C3 can also be negative. In this case, the cubic has one real and two 
imaginary roots. Defining 

1/2 ^ \ 1/3 



A 




we can write the real root as 



z = A- 



■i\cM 




{Cz < 0). 



(A2) 
(A3) 
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